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l0,a ." y ;, nle l r l nal,y J '? n « tin « "ncentrators (DTIRCs) which can achieve concentrations close to the 
thermodynam.callv allowed l.m.u are introduced. Ceneral design methodologies are given explicitly Geo 
• ^metrical and optical properties of DTIRCs are discussed. explicitly. Oeo- 



I. Introduction 

A. General 

An optical concentrating system can be character- 
ized by its angular field of view and its entrance and 
exit apertures. The geometric concentration defined 
as the area ratio of entrance to exit is limited by the 
angular field of view. In a 3-D geometry, the maxi- 
mum concentration (thermodynamic limit) for an op- 
tical system with a conical angular field of view of 
(half-acceptance angle) ±6 a is given by 



n;/n|sin 5 « 0 . 



(1) 



where n, and n 2 are the refractive indices of the media 
in which the entrance and exit apertures are immersed 
respectively.' It is also known that the concentration 
of most simple conventional optical systems (e.g lens 
mirrors) are far below these limits. An //stop Fresnei 
lens falls short of these limits by a factor of Af 7 in 
vacuo. 2 

In the last two decades, intensive research has been 
carried out in designing concentrators (in both 2-D and 
3-D geometry) which have concentrations achieving or 
approaching these limits. 1 - 3 A well-known example of 
such a concentrator is t he compound parabolic concen- 
trator (C PC) which has been used in solar energy con- 
c entration systems 4 and high-e neit-v Clieienkov l. P h r ■ 
detect" - " ' " — — 



In thjs paper, we will discuss a new family of optical 
elements, dielectric totally internally reflecting con- 
centrators (DTIRCs). DTIRCs combine the front 
surface refraction with total internal reflection from 
the sidewall to achieve concentrations close to the 
theoretical maximum limits. Compared with CPCs 
because of the larger refractive index of dielectric and 
the use of a curved front surface, DTIRCs have two 
advantages: higher concentrations and smaller sizes, 
iney can be used in many areas, e.g., solar energy 
concentrators and fiber optics, where the concentra- 
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tion or ef ficjepl transport o£ light is required, ror 
example, the dielectric CPC (DCPC) h^sbein~devel- 



For 



oped as a single-element photovoltaic concentrator. 6 - 7 
In Sec. II, we will discuss the general design method- 
ologies for DTIRCs. Their geometrical and optical 
properties will be analyzed in detail and comparisons 
to experimental data made in the following sections. 

B. Nomenclature 

6 0 — acceptance angle; 
d\ = entrance aperture; 
do = exit aperture; 

<fi = front surface arc angle; 

n = index of refraction. 

II. General Design Method of DTIRCs 

When light rays are incident on a boundary of a 
medium of higher refractive index from one with a 
lower index at an angle greater than the so-called criti- 
cal angle, the rays will all be reflected back to the 
region with a higher index. This physical phenome- 
non is known as total internal reflection (TIR), and the 
critical incident angle B c is derived from Snell's law: 



e t = sinfn./n,). 



(2) 



300 



APPL.EDOPT.es / Vol. 26. No 2 / ,5 January ,987 



where n, and n ? are refractive indices and n, < n->. 

It is possible to design a solid dielectric concentrator 
in which the condition for TIR is fulfilled at the side 
profile boundary for all rays incident on the concentra- 
tor within the design acceptance angle. Since it uses 
total internal reflection properties of the dielectric 
medium, it is known as DTIRC, A DTIRC consists of 
three parts: the front surface which may be curved- 
the totally internally reflecting sidewall (profile)- and 
the exit aperture. Only rays within the design accep- 
tance angle (corresponding to an angular field of view) 
need reach the exit aperture. Most rays incident on 
the front surface beyond the acceptance angle will 
eventually exit from the side profile. Ravs incident at 
the acceptance angle are known as extreme rays, and 
they play an important role in determining the side 
profile of DTIRC. As for CPCs, we will solve the ->-D 
case and form 3-D DTIRCs by rotating the 2-D profile 
about its symmetry axis. In this paper, we will princi- 





: :.L S.dev.ewofaDTIRC. The extreme rays hitting the por- 
-aPuPi converge into/' ,. Rays hitting the portion P 2 J> 3 can exit 
.« -vnatever angle barely satisfies TIR (the maximum concentration 
method) or in parallel (the phase-conserving method) 



pally be concerned with the properties of the 3-D solu- 
tions. 

Figure 1 shows an example of a DTIRC A set of 
r -ys incident on the front surface at the extreme angle 
it shown. Each passes through the front curved sur- 
l.ice, is refracted, and subsequently is totally internally 
:m lected from the sidewaJl and eventually hits the exit 
aperture. In determining the side profile, we maxi- 
•mze the slope of the side profile curve consistent with 
rejecting the extreme entrance rays after refraction by 
the curved front surface onto the exit subject to various 
Mitisidiary conditions we may wish to impose For 
example in addition to TIR, we may require rays to 

6Ml W nii n n o angle n0t exceedin & a predetermined an- 
Sie. DI IRCs generated in this subsidiary condition 
-an be used m fiber optics where rays must enter fiber 
at angles smaller than a certain critical angle for trans- 
mission. 

We will discuss two design methods for DTIRCs. In 
jne first method, referred to as the maximum concen- 
tration method, there are no subsidiary conditions oth- 
Zu J R ,m P° sed . ar >d it yields the maximum 
achievable concentration with TIR. In the second 
method, the exiting extreme rays are parallel, that is 

ney form a new wavefront. This approach is referred 
'o as the phase-conserving method. 

Let us consider the example shown in Fig. 1. For 
-implicity (in both designing and manufacturing), we 

"i restrict ourselves to the case where the front sur- 
ace ,s a portlon of a sphere , n gcneral> wg Cfln divjde 

ne side profile into two parts: the portion from P, to 
and the portion from P 2 to P 3 . Between P, and P, 

sin-le Tm "i^"* di [ eCted t0 the COrner P 3 after a 
Tin ! a T he ra ? hittin B P 2 J' u9t bare| y satisfies 
"< and exits from P,. Between P 2 and P 3 , the refrac- 

extrem ,S ^ h ' eh en ° Ueh 10 totally reflect the 
profilp in?!? 5 v There are two Wfl y s to desi 'g n the 
be mlin? t . h,s J re * ,on - One can simply require that TIR 
whIE * T d 3nd allow the Meeting rays to assume 
exit ,n , l T reSult9 as ,on S as th<f y strike the 
an-les f„ u " ThlS results in a distribution of exit 
- lor the extreme incident angles and develops 




hit 

Aptrtui* 



-rl»»l W*v*(tont 

Fig. 2. Typical optical path conaisla of four parti: /,,/,, I, and /« 
tor rays converging at the bottom left-hand corner, the portion /< i, 

zero. 



the maximum possible concentration achievable with 
an edge ray solution incorporating TIR. The second 
method is to slope the profile so that the reflecting rays 
are all parallel to one another with a direction corre- 
sponding to the line from P 2 to P 3 . This means that 
the extreme rays reflected from the portion between P 2 
and P 3 emerge with a well-defined wavefront 

The two methods yield solutions which are very 
similar geometrically, but the advantage of the latter is 
that an analytical solution can be determined That 
is since extreme rays reflected from portion P { J> 2 of 
the profile converge onto a point and rays from the 
portion P 2t P 3 exit in parallel. Fermafs principle can 
be used to generate the coordinates of the side profile 
In general, the total optical path length of a ray from 
the entrance wavefront to the exit wavefront consists 
of four parts: /,, l 2 , / 3 , and /, as shown in Fig. 2. The 
initial wavefront is shown just tangent to the curved 
front surface. The portion /, is from a point on the 
initial wavefront to the surface of DTIRC, /, is from 
there to the reflection point, l 3 is from the reflection 

Km? * /I' ? °f DTIRQ 3nd U is from the exit of 
P / is zero Wavefront - For ra >'s converging to 

According to Fermafs principle the total optical 
path ength is a constant for every ray connecting the 
initial and final wavefront. Thus we have 

/ nd* =/, + „(/, + /, + /,, = nmsl . , 3) 

Combining Eq. (3) with DTIRC height and length 
constraints, we can solve for /,, /,, an d /«. Then the 
profile coordinates can be calculated analytically. 
1 he detailed relationships are derived in the Appen- 
dix. With the maximum concentration method, rays 
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from portion P X JP 2 are still convergent to P 3 . Thus 
the method outlined above is still applicable. But for 
portion conserving the optical path can no long- 

er be applied, since there is no well-defined wavefront. 
Recursion relations are developed to calculate the co- 
ordinates of the side profile. The detailed formula- 
tions are also shown in the Appendix. 

We have developed computer programs to generate 
the side profile curves numerically. The programs 
require as inputs the front surface arc angle, the refrac- 
tive index of the dielectric, the exit aperture, and the 
acceptance angle. The programs assume a trial en- 
trance aperture and, therefore, trial height, then step 
through a set of extreme rays (where the total number 
of extreme rays depends on the precision requirement) 
to calculate coordinates of the DTIRC profile. It be- 
. gins with the extreme ray reflected at P 3 and ends with 
the extreme ray entering the edge of the entrance 
aperture at P,. The program then compares the en- 
trance aperture generated with the trial aperture. 
The difference between two entrance apertures is used 
to adjust the trial aperture and height, and the proce- 
dure iterates until the trial aperture and calculated 
aperture converge to within a predetermined error 
tolerance. The entrance aperture is then defined, and 
X-Y coordinates of the DTIRC are calculated. The 
number of extreme rays determines the number of 
coordinates and also the numerical precisions. The 
following discussions are based on numerical solutions 
for DTIRCs which have 200 pairs of X-Y coordinates 
for the side profile. 

When the radius of curvature of the spherical sur- 
face is small, we can treat the surface as an ideal lens 
which images the extreme pencil onto a point. In this 
approximation, for DTIRs designed using the maxi- 
mum concentration method, the portion between P u 
P2 is a hyperbola, while the portion between P 2r P 3 is an 
arc of an equiangular spiral (see Fig. 3). 

III. Geometrical Properties of DTIRCs 

As indicated above there are four degrees of freedom 
in designing a DTIRC. Among these, the refractive 
index and the exit aperture are fixed in most situa- 
tions. Thus we will focus our discussion on the effects 
caused by variations on the acceptance angle and front 
arc angle. It is found that the geometrical concentra- 
tion defined as the area ratio of the entrance to exit 
depends not only on the acceptance angle but also the 
front arc angle, and this holds for DTIRCs designed 
using either method. In the phase-conserving meth- 
od, we require not only that extreme rays satisfy TIR 
but also exit in parallel. Thus with the same front arc 
angle, DTIRCs generated using the phase-conserving 
method have a slightly lower concentration than 
DTIRCs using the maximum concentration method. 
Table I compares the concentrations for DTIRCs gen- 
erated using the two methods at various acceptance 
angles. 

The TIR condition is more difficult to satisfy as the 
front surface becomes more curved because rays are 
bent more sharply. Thus the concentration becomes 




Fig. 3. When the front surface curvature is small, (he portion P tw P 7 
is just a hyperbola, and the portion P 7 ,P 2 is JU st an arc of an equian- 
gular spiral. 



Table I. Companion of Concentration* (Index = 1.5, A/c Anglo « 30°. 
_____ E»H Aperture = 1 cm) 



Acceptance 
angle 

(deg) 


Phase- 
conserving 
method 


Maximum 
concentration 
method 


18.0 


21.81 


22.09 


20.0 


17.56 


17.81 


22.0 


14.44 


14.67 


24.0 


12.11 


12.32 


26.0 


10.24 


10.43 


28.0 


8.76 


9.00 



Table II. Concentration* va Arc Angle <lndei * 1.5, Acceptance Angle - 
22°, Eilt = Um) 



Arc 
angle 
(deg) 


Phase- 
conserving 
method 


Maximum 
concentration 
method 


10 


15.37 


15.44 


20 


14.98 


15.05 


30 


14.44 


14.67 


40 


13.84 


14.36 


45 


13.47 


14.29 


50 


13.10 


14.21 



smaller for a DTIRC as the front arc angle is increased. 
Table II shows the concentrations for DTIRCs with 
different front surface curvatures for a fixed accep- 
tance angle for both maximum concentration and 
phase-conserving methods. From Tables I and II, it is 
obvious that DTIRCs generated using the maximum 
concentration method have slightly higher concentra- 
tions. DTIRCs designed using this method are pre- 
ferred in applications where concentration is the only 
concern, i.e., solar energy photovoltaic concentra- 
tors.- 3 Since for DTIRCs designed using the phase- 
conserving method, the maximum exiting angle is ex- 
pressed in terms of the four input parameters 
analytically, we can choose appropriate parameters in 
designing DTIRCs so that no ray will exit with an angle 
larger than a predetermined angle. DTIRCs with a 
maximum predetermined exiling angle can be used in 
fiber optics where the angular requirement is crucial. 
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t -<z. 4. Concentrations as functions of one acceptance angle for 
V.ih DTIRC and CPC. With the same acceptance angle, DTIRC 
has a higher concentration. 




ACCEPTANCE ANQLE fOfCUtt) 

•=>. Idealness of the DTIRC decreases with the acceptance 
angle, while the idealness of CPC is one. 



^pite all these comments, the geometrical differ- 
ences between DTIRCs designed using different meth- 
ods are quite small. Especially for the DTIRCs with a 
-mailer acceptance angle and/or smaller front curva- 
tures, the difference is insignificant in practical appli- 
cations. In the following discussions, we will only 
ronwder DTRICs designed using a maximum concen- 
tration method. 

Figure 4 shows the geometrical concentrations of a 
mirror-constructed CPC and DTIRCs designed for 
various angular acceptances. If one defines idealness 
s me ratio of geometrical concentration of a DTIRC to 
ne maximum limit given by Eq. (1) for a given accep- 
ance angle, a quantitative relationship between ideal- 
ness and the acceptance angle for a DTIRC is shown in 
for PPr- • lnce J t0lal internal reflection is not required 
tha, ,k ' US ' dealness is or, e. We can see from Fig. 5 
DTIRP *?* £' the acce P tance angle is, the closer the 
th» 'v° tne thermodynamic limit. Unlike CPCs. 
^ne curved front surface of DTIRC functions like a lens 
n.ch concentrates before the rays are reflected. Ap- 

*ho te t the D-FlRP ^ T ' ^1^' the 

DTWr ine L D , IIRC - F'gure 6 shows a family of 

given ^rent front surface curvatures for a 

with * n°? P , tanCe an - le " A ^electric CPC (DTIRC 
surf, ° nt entran ee) is also shown. The front 

2Sff iR ? ar vr h smal,er - it ™ d "- 

the xnn lr.u 9 he ' ght as lhe dlsta "ee from the exit to 

to the ?n V u C Pf0file flnd total hei S ht from l he exit 
l °P ol the curved front surface, quantitative 




Fig. 6. Family of DTIRCs with a different front surface arc angle 
Because of curved front surfaces. DTIRCs can be made smaller. 
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Fig. 7. DTIRC height and the total height'aecrease with the front 
surface arc angle. 

relations between heights of DTIRCs and front arc 
angles are shown in Fig. 7. 

We should point out that as the front surface varies 
%°S£ln I curvat f ure t0 a ,a 'ger one, the side profile 
oi u l IKC changes from convex to concave. This sug- 
gests the possibility that with a certain front curvature 
the side profile of DTIRC can be well approximated by 
a straight line In fact this turns out to be the case. 
Straight side DTIRCs are expected to be easier to 
produce Table III shows the linear correlation coeffi- 
cient r of the side profile of DTIRC at various front 
curvatures. 9 As one can see from this table, the side 
profile of DTIRC can be very well approximated by a 
straightlineforafrontsurfacearcangleof~50° Ifwe 
relax the requirement that the front surface he spheri- 
cal, it is possible that the side profile is exactlv a 
straight line in certain circumstances. Aspherical 
front curved DTIRC will be discussed in a later paper. 
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T«W. Ul. Arc An 8 ,. „. LlnMr C or,.,.„ 0 „ Co.n,c..n. (l nd.« ol R.„ac,.on 
- - 1.47, Acceptance Angle « 22°) 



Arc angle (deg) 



Coefficicnl r 



10 
20 
30 
40 
45 
50 
55 
60 
65 



0.95910 
0.98268 
0.99326 
0.D98U 
0.99925 
0.99966 
0.99945 
0.99862 
0.99723 
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axgle (Otonet) 

Comparisons of compu.er simulated angular acceptance 
'.lh data for a DTIRC with a design acceptance angle of 
2\.S4'. 



IV. Optical Properties of DTIRC 

The optica] performance of DTIRC can be best de- 
scribed by its angular acceptance. The angular accep- 
tance curve IS defined as the normalized optical 
throughput vs light incident angle 8. The angu a 
acceptance curve for DTIRC is very similar to Cp£ 
and has a clear cutoff at the designed acceptance angle 
i-'gure 8 shows a computer simulated angular accep- 
tance curve for a physically ideal DTIRC designed for 
an acceptance angle equal to 21.64°. The roundine 
and slopmg falloff of the ideal DTIRC in the m £ uEJ 
acceptance curve is due to skew ray losses. The angu- 
lar acceptance can also be determined experimentally 
We measured a DTIRC designed with an acceptance 
angle of 21.64-. front arc angle of 65% and an ex" 
aper ure of 1 cm. A 1-cm diam solar cell optically 
coupled to the exit of the DTIRC was used as the 
dZt th Tt hpUt detector - The experiment was 

drcui curtt' ° n 3 dear d3y - BeC3USe the ce " sho " 
circuit current is proportional to the amount of light it 

receives, by measuring the cell short circuit current a 

different angles, the angular acceptance can be deter 

mined a ter the effects of diffuse sunlight are con c - 

incident an.T/ "t lhrOU S h ? ul F » ■ function of 
incident angle 0 can be expressed as 

mal/zed' to Sh ° rt C , irCUil CUfrent with DTIRC nor- 
drcu itcurren? ms ° ,at ' on - l > the cell alone short 
circuit current normalized to total insolation, and C g is 



lieTnThScat. 00 """^ 31 ' 011 ° f ^ ° TIRC ' which i£ 

la4 C t 0 h7eff ter , m0 f d r el ^ 3,50 devel «P^ which calcu- 
n thJ A ! , CtS ° f f !, 0nt Surface ^flection, absorption 
in he dielectric ;. and the reflections at the dielectric- 
cell interface. Figure 8 also shows the comparisons of 
computer simulations with the measured da'ta From 
this figure. ,t , s dear that theory and experiment agS 

V. Conclusions 

con\en e tT fl t C or S H S °- l° ta] internal election 

concentrator designed according to optical principles 

ha concentrations approaching the theoretical rZ* 

mum hm.U. Compared with CPCs. DTIRCs offer 

higher concentrations and smaller physical sizes 

which a re .mportant in practical applications. Since 

SSTK COnd '^ Can be ira P° sed on designing 
DTIRCs. they can be designed to meet various kinds of 
requ.rements m situations where concentrations are 
necessary. Applications of DTIRCs as a second^ 
£S? U r°e r m Ph ° t0v0ltaic wi 'J ^ discussed 

This work is supported in part by the U.S. Depart- 

?z:tz' m developme part ° f ihe —k—tai 

Appendix: Mathematical Formulations 

A: Phase-Conserving Method 

Given a wavefront of extreme rays, the part of it 
which hits portion P } ,P 2 of the side profile of DTIRC is 

3£?ft lnt ° b D 0t o° m C ° rner P 3- Another par 
which hits portion P 2f /> 3 must leave the aperture as a 
new wavefront. If we use /, (f = 1,2.3,4) for different 
par s of the optical path shown in Fig. 2. then. accorS- 

+ Vu I'?! 1 ,* P" nc,p J e ' tne optical Path length C = /, 
+ n(/ 2 + l 3 + / 4 ) must be a constant 

The radius of front curvature R can be expressed in 
terms of front arc angle * and entrance diameter d by 
the following equation: . 

ft = d,/(2siivl. (A1) 

Now consider an extreme ray which hits the front 
surface corresponding to an arc angle 0 shown in Fig 2 
This ray will leave the front surface with an angle 8' 
relative to vertical axis: g 

= sin-'|si,i«l + « o )/ n | -h. 

And the optical path /, is given by 

/, = 2H>in : [ui + 

DTiRCh 6 f K C , l w hat ?V y fr ° m Pl ?hould hit 
DTIRC height H can be determined. Using Eq. (A2 

it can be shown that 

H = + d j ,.,„„„_ (A4) 

where 0 = maxW) = R in-|sin(fl 0 - ,,)/„) + c , and the 
optical path C is given by ' ™ 



(A2I 



(A3) 
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C = 2R s\n 2 \(8 0 + v,)/2] + n(d l + <f 0 }/[2 sin(O)]. (A5) 

If 6 is such that the ray hits portion P X JP 2 of the side 
profile, the ray goes into P 3 and / 4 = 0. We have three 
unknowns (kJ^hy) and three equations: 

Optica) path length: C - / 1 + n\l 2 + (/^ + l\ y )^\- ( A€) 
Height: / 2 cos(*') + l Zy - R[cos($) - cosM] - H\ (A7) 
Length: R sin(0) + / 2 sin(0') - ~ -djl. (A8) 

From these equations l 2 can be solved: 

U = |[(c - l x )/n} 2 - [R sin(€) ^d^] 2 - H 7 ]/ 

\{c - l x )fn + [R ain(0) + dJ2] sinfl' - // cos*'|. (A9) 

The extreme rays hitting a portion of the P 2r P 3 exit 
with the same angle 6 Qt which is determined by requir- 
ing that the ray from Pi to P 3, just barely satisfies the 
TIR condition at P 2 . And the exit angle 8 0 is given by 
the following formula: 



* 0 - r - e - 2* 



(A 10) 



Where 6c is a TIR critical angle given by Eq. (2). 

Once again we have three unknowns (l 2y l 3 J 4 ) and 
three equations: 



Optical path length: C - /, + n{l 7 + / 3 + /<); 
height: / 2 cos(0') + / 3 cos<* 0 ) « //; 
U-ngth: /? sin(0) + l 7 sin(0') - / 3 sin(0 o ) 



(All) 
(A12) 



- - V2 + l A /sin(B 0 ). (A 13) 

To solve these equations, we obtain 

- \(C - /,)/2 - H/cos(6 0 ) - (/? sin(0) 
- H/cos(6 0 ) + sin(0 o )|/ 

[1 - cos(0')/cos(0 o ) + sin(tf') sin{^ 

+ cos(0') sin(0 o ) tan(0 o )], (AM) 

'3 " |rY - / 2 cos(0'))/cos(0 o ), (A15) 

- I* sin(fl) + / 2 sin(ff') - / 3 sin<0 o ) + d<j2\ sin(* 0 ). (A16) 

Once we have / 2 , then x and y coordinates are given as 
follows: 

x = R sin{0) + / 2 sin(0'); (A17) 
y = H-l 2 cos(6') + /?[cos(0) - cos(*)J. (A18) 

B: Maximum Concentration Method 

Compared with the phase-conserving method, the 
Maximum concentration method does not require rays 
hitting a portion P 2 ,P 3 of the side profile exiting with 
the same angle. For rays hitting a portion P lt P 2 , the 
methodology is the same for both phase-conserving 
and maximum concentration. Thus the formula de- 
rived above for the portion P lf P 2 [Eqs. (A1)-(A9)] is 
st »N applicable. We will focus on discussions on the 
Portion P 2 J>. d where there is no longer a clearly defined 
w av e front. 



Assume that a point with coordinates (X lt Y t ) in the 
portion P 2 JP 3 has been determined. .The next point 
(X^uYt + i) is just the intersection of the next extreme 
ray with a straight line extending from (X„Y,) with 
such a slope that this ray will just be totally internally 
reflected. It can be shown that (X^ U Y^\) is connect- 
ed to (X it Y t ) by the following recursion relations: 

+ , = + sin(0) - V, + X, tan(0 f + $') - R s j n (0) cot(0')]/ 
\l&n(8 c + $') - cot(0')]; (Bi) 
y t + i = (X,^- X t ) tantf, + 8') + V,. (B2) 

Let us use TV as the number of total extreme rays which 
determines the solution precisions. In the computer 
program, one starts with P 3 , which has coordinates 
(0,d 0 /2). The next point is determined by considering* 
the next extreme ray, which hits the front surface 
corresponding to an arc angle 6 = -(N - l)//v>. Then 
the next extreme ray is traced in the same fashion unt* 1 * 
a ray exits P 3 . Then the phase-conserving formula 
tions are used. 
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